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Graph-Directed Fractal Interpolation Fnnctions 


Ali DENiZ * and Yunus OZDEMIr t 


Abstract 

It is known that there exists a function interpolating a given data set such that 
the graph of the function is the attractor of an iterated function system which is 
called fractal interpolation function. We generalize the notion of fractal interpolation 
function to the graph-directed case and prove that for a finite number of data sets there 
exist interpolation functions each of which interpolates corresponding data set in 
such that the graphs of the interpolation functions are attractors of a graph-directed 
iterated function system. 
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iterated function system 
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1 Introduction 

Barnsley introduced fractal interpolation functions (FIF) using iterated function system 
(IFS) theory which is an important part of fractals (see [I], [2], [11] for further information). 
He showed that there exists a function interpolating a given data set such that the graph 
of this function is the attractor of an IFS. In recent decades it has been widely used in 
various fields such as approximation theory, image compression, modeling of signals and 
in many other scientific areas (a,Ha,[10],ID). 

Let us summarize the notions of IFS and FIF firstly. An IFS is a finite collection 
of contraction mappings ipi : X ^ X (i = 1, ...,n) on a complete metric space X. It is 
known that there exists a unique nonempty compact set A satisfying 

n 

A = \J ipi(A) C A 

i=l 

‘Department of Mathematics, Anadolu University, Eskisehir, Turkey. Email: adeniz@anadolu.edu.tr 
^Department of Mathematics, Anadolu University, Eskisehir, Turkey. Email: yunuso@anadolu.edu.tr 


1 



which is called the attractor of the IFS (see i)- 
A data set is a set of points 

2) = {(xo,-Fo), (a;i,-Fi),..., {xn,Fn)} C 

such that xo < < • • • < xat with N > 2. An interpolation function corresponding to 

the data set 2) is a continuous function / : [xq, xat] —>■ M whose graph passes through the 
points of 2), i.e. /(xj) = Fi for alH = 0,1,..., iV. 

Barnsley constructed an iterated function system whose attractor is an interpolation 
function corresponding to the data set 2) as follows: 

Let cji : —7> be affine functions of the form 



where the Oj, Ci,di, e*, fi are real numbers for i = 1,2,..., N. Using the constraints 

uJi{xo,Fo) = {xi-i,Fi_i) and uJi{xN, Fn) = {xi, Fi), 

one can obtain the coefficients ai,Ci,di,ei, fi easily in terms of data points {{xi, Fi)}^^ 
and the parameters di (so-called vertical scaling factor) for i = 1,2, ...,N as 

Xi Xi—\ 

Oj = - 

XN — Xo 

X NXi—\ XQXi 

Gi = - 

XN - 

_ Fi — Fj_i F/v — -^0 
Ci — di 

XN — Xo XN — Xo 

XNFi_i - XoFi XnFq - XqFn 
Ji — di 

XN — Xq Xn — Xq 

(see m, m for details). 

In [21 Teorem 2.1, p.217], it is shown that cui,..., cjat} is an iterated function system 
(associated with the data set 2)) in case of \di\ < 1 for z = 1, 

The attractor of the IFS ...,a; 7 v} constructed above is the graph of an interpo¬ 

lation function which is called (affine) fractal interpolation function (FIF) corresponding 
to the data set 2). Note that there exist infinitely many fractal interpolation functions 
depending on the parameters di with |dj| < 1 for i = 1,..., A^. 

Example 1. Consider the IFS (constructed as above) assoeiated with the data set 

{(0,0), (3, 5), (6,4), (10,1)} 
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where we choose di = 0.25, (i 2 = 0.5, da = 0.25. The IPS consists of the contractions 

u!i{x,y) = (0.3x, 0.475x + 0.25y) 
uj2(x,y) = (0.3x + 3, —0.15x + 0.5?/+ 5) 
uj3{x,y) = (0.4x + 6, —0.325x + 0.25?/+ 4) 

and its attractor is a FIF whose graph depicted in Figure [3 



Figure 1: The attactor of the IFS cji,a; 2 , ws} given in Example [3 


There is a generalization of iterated function systems, the so-called graph-directed 
iterated function systems (GIFS). In this work, we generalize the notion of fractal interpo¬ 
lation function to the graph-directed case. In modelling with fractal interpolation, graph 
directed fractals can have some advantages since they give us a chance to code several 
data sets in an interpolation function. An interpolant obtained as graph directed fractal 
interpolation function contains some information from each of the data sets under suitable 
conditions. 

We briefly summarize the graph-directed iterated function systems. 

Let G = {V-, E) be a directed graph where V is the set of vertices and E is the set of 
edges. For a,j3^V, the set of edges from a to /3 is denoted by E^h and its elements by 
i = 1,2,..., where iF"^ denotes the number of elements of E°‘^. We assume that 

E" := IJ ^ 0. 

p&v 
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Let {X°‘ I a G y} be a finite collection of complete metric spaces and —>■ X°‘ 

be a contraction mapping corresponding to the edge (but in the opposite direction of 
e“^) for i = 1, 2,... , K°‘^ and a, /3 G IL. It can be shown that there exist a unique family 
of nonempty compact sets C X" such that 

= U U 

P&v i=l 

(see IS]). The system |x"; is called a graph-directed iterated function system (GIFS) 
realizing the graph Q and A" ’s (a G V) are called the attractors of the system (see [5],[8| 
for more details). 

Example 2. Let Q = {V,E) be a directed graph with V = {1,2} as shown in Figure\^a). 
Let X^ be the unite square and X‘^ be the right triangle as indicated in Figure lE(b)- 
Consider the contractions —)• X" (a,f3 ^ V, i = corresponding to 

each edge of the graph (in the opposite direction). A pictorial description of the maps are 
shown in Figure\^c), the explicit expressions of them can be found in The attractors 
of this graph-directed system realizing the graph Q look like as given in Figure\^d). 


2 Construction of Graph-Directed Fractal Interpolation Func¬ 
tions 


In this section, we will construct a graph-directed iterated function system for a finite 
number of data sets such that each attractor of the system is the graph of an interpolation 
function (which interpolates a corresponding data set) what we call a graph-directed fractal 
interpolation function. 

Let 

= {(4, Fo“), (x?, Ff),..., (x^^, F^Jj (1) 


be data sets in with > 2 for all a = 1,2,... 
satisfy 

-A -L_L < 1 

-^0 

for all a 7^ /3, a, /3 = 1,..., n and j = 1, ..., Np. 


n. We assume that these data sets 

( 2 ) 


Remark 1. The condition on the data sets will provide that the maps to be constructed 
further in this section are contractions. Note that if the original data sets do not satisfy 
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Figure 2: (a) A directed graph G = {V,E) where V = {1,2}. (b) The spaces and X‘^. 
(c) Pictorial description of the contractions ojf’^ which is determined by a copy of X^ in 
A" for all a,j5 i = 1, . (d) Sixth stage of the iteration of the GIFS. 


the condition one can make the data sets to satisfy the condition by adding new points 
artificially. Instead of changing the original data sets by adding artificial points, we prefer 
the data sets to satisfy the condition 

Now we can state the following theorem as the main result of this paper. 

Theorem 1. Let D" = {{xq,Fq), ..., be data sets in for a = 

satisfying Then there exists a GIFS, with attractors A" (a = 1,2,... ,n), 
such that A“ is the graph of a function which interpolates the data set ID" for each a. 

We call an interpolation function whose graph is the attractor of a graph-directed 
iterated function system a graph-directed fractal interpolation function. 

In the rest of this section we deal with the construction of a GIFS satisfying Theorem [1] 
and for simplicity we consider the case n = 2. In Section [3l we prove that Theorem [1] 
holds for a GIFS to be constructed in this section. 
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Let 


{(xo, -Fo), (xi, Fi), ..., (xat, Fn)} 

{{yo,Go), (?/i,Gi),..., (?/m,G'm)} 


2)^ = 
1)2 = 


be two data sets in with N,M > 2. In the case of n = 2 the condition ([2|) reduces 


and 


Um - Vo 


XN — Xo 


( 3 ) 


for all z = 1,2,..., iV and j = 1, 2,... , M. 

Clearly, one can form a directed graph Q = {V, E) with V = {1, 2} such that 

+ = N and = M (4) 


(recall that is the number of elements of for a,/3 G V). Such a graph is shown 
in Figure [3l 



Figure 3: A directed graph {V, E) where V = {1,2}. 


Remark 2. For the general case, the eondition & on the number of edges between the 
vertiees of the graph will be 

A"! + + ... + A"” = Na 

for a = 1,2,... ,n. 

Remark 3. Note that we allow the ease A"^ = 0. This may be lead to a situation that 
a data set ean not he coded in an interpolant. If one eonstruet a strongly eonneeted graph 
(i.e. there is a directed path from the initial vertex a to the terminal vertex (5 for each 
pair a,l3&V of the graph) satisfying the eondition 0 then each interpolant will nestle 
some information from each of the data sets and thus the graph-directed self-similarity of 
the interpolants become more genuine. 
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To obtain a GIFS associated with the data sets 0“, (a = 1, 2) and realizing the graph 
Q consider the afhne functions —>■ i = 1, 2,... , of the form 


03 , . I “f 0 



where the a“^, c“^, e“^, are all real numbers for a, (3 £ {1,2}. 

We constrain these functions as follows: 


wI^{xn,Fn) = {xi,Fi) 


for i = 1,2,..., F 


11 


WilKiiiyo^Go) = {xi-i,Fi_i) 
= {Xi,Fi) 

w^^{xo,Fo) = {yi_i,Gi-i) 


for i = + 1,..., +K^'^ = N 


11 


'12 


wf-{xN,FN) = {yi,Gi) 


for i = 1, 2,... , F 


21 


WilK^iiyo,Go) = {yi-i,Gi-i) 
wfij^2i{yM,GM) = {yi,Gi) 


for i = + 1,..., + K'^'^ = M 


21 


'22 


Note that each constraint should be applied in the case of K°‘^ > 1. 

Remark 4. Complying with the mapping relationship determined by the graph Q, one can 
rearrange the maps in an arbitrary order. Since our aim is to show the existence of a 
graph-directed interpolation function, we arrange the mappings according to the order of 
data sets. 

From each group of above conditions we get the following equation systems respectively: 


Xi-i = al^xo + e, 

F,_i = cj^xo + dj^Fo + fl^ 

Xi = a\^XN + 

Fi = cj^XN + dj^FN + fl^ 


For all i = 1,2,..., F 


11 


(5) 


For alH = F^^ + l,...,F { 


Xi-l 

F^-l 

Xi 

Fi 


„12 ,, _i_ „12 

c]^K^iyo + difj^iiGo + 

allK^iVM + 

c]‘^K^iyM + dffj^iiGM + fl^K^i 


( 6 ) 
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Vi-i 


For alH = 1,2,..., 


21 


= af^xo + 


Gi-i = 

Vi = 
= 


I 


i 

' 21 ; 


d^xo + di^F^ + fl^ 
af-XN + ef^ 
cfxN + dfFN + /f 


(7) 


For all i = i^2i + 




1, — /,22 ,, I „22 

Vi-l — O'i_i^2iyo + 


G,_i = + df_^2iGo + 

Vi = afip.2iyM + edj^2i 

Gi = cJ^j^2iyM + ddj^iiGM + 


( 8 ) 


Choosing as a parameter, one can easily solve the linear equation systems (|^~(|5|) 
to obtain a“^, c“^, e“^, /“^ for a, /3 € {1, 2}, i = 1, 2,... , in terms of the data points 
and df^ as follows: 


^i—1 

af^ = - 




4^ = 


XN — Xo 
X NXi—\ XqXi 
XN — Xq 
Fi — Fj i 


-dd 


XN — Xo 


Fn — Fq 

XN — Xo 


12 ^i—1 

ar = - 


yM - yo 

12 yMXi—\ 

ef =- 


yoxi 


cF = 


yM - yo 

Fi -Fi—1 ^12^^ ^0 

yM - yo 


VM - yo 


yii_ xnFi—x xoFi ^^-^xnFo xq-F/v j,x2 _ yMFi—x yoFi yoGM 


XN — Xo 


XN — Xo 


VM - yo 


VM - yo 


ar = 


Vi yi—1 
XN — Xo 


xoyi 


21 XNyi—1 

= - 

XN - Xo 

^21 _ Gi Gi—X ^21 Fn Fq 


XN — Xo 


XN — Xo 


j 2 X _ XNGi—X XoGi ^21 XnFo XoFn 


XN — a:o 


XN — Xo 


„22 _ Hi yi-1 


yM - yo 
^22 _ yMyi—1 — 


yoyi 


yM - yo 


^22 _ Gi—x ^ 22 ^^ ^0 


yM - yo 


yM - yo 


f 22 _ yMGi-X — yoGi j 22 yMGo — VoGm 

h — ~ 


VM - yo 


yM - yo 


To show that the obtained system is a GIFS, we need to ensure that is 

a contraction for each i, a and /3. Although ’s need not be contractions with respect 

to the Euclidean metric on M^, the following lemma guarantees that they are contractions 
with respect to a metric which is equivalent to the Euclidean one. 




























Lemma 1. Let he the ajfine functions constructed above associated with the data sets 
S", (a = 1, 2) which satisfy Let \df^\ < 1 for all a, (3 G {1, 2} and i = 1,2,..., 

Then the system is a graph-directed iterated function system associated with 

the data sets D°‘, {a = 1,2) and realizing the graph Q. 

Proof. Condition ([3|) and [21 Theorem 2.1] make each of a contraction with respect 
to a metric which is equivalent to the Euclidean metric. See jH Theorem 2.1] for more 
details. □ 

Depending on the parameters we have infinitely many GIFS realizing 

the graph Q. 

Example 3. Consider the GIFS associated with data sets 

= {(0,5), (1,4), (2,1), (3,1), (4, 4), (5,5)} 

= {(0,1), (1,2), (2, 3), (3, 2), (4,1)} 

realizing the graph with = 3, = 2, = 1, = 3 as shown in Figure [2[ 

We choose dff^ = - for all a, ft G {1,2} and i = The attractors and A'^ 

O 

of the system are given in Figure A^ (resp A"^) is the the graph of a function which 
interpolates (resp. D'^). 

We also present the attractors of another GIFS which is obtained by changing only the 
vertical scaling factors as = d^^ = d}^^ = 0.25,= d^^ = \,d^ = \,d‘^ = d'^ = 
dg^ = ^ (see Figure\^. 



Figure 4: The directed graph used in Example [3l 
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Figure 5: The attractors (left) and (right) of the first GIFS given in Example [3l 




Fignre 6: The attractors of the second GIFS given in Example [3] with different vertical 
scaling factors. 
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3 Proof of the Main Theorem 


In this section, we prove that a GIFS constructed in Section [2] satisfies Theorem m As 
mentioned before, we consider the case n = 2. 

Specifying the coefficients such that < 1 for all a,/I € 1/ = {1,2} and 

i = 1,2,... let be a graph-directed iterated function system associated 

with the data sets S", {a = 1,2) constructed in Section [2j We prove that its attractor 
system {A^,A‘^} is a pair of compact subsets of such that A^ (resp. A^) is the graph 
of a function which interpolates the data set (resp. A). 

Let 

•^ = {/ I / : [3;o, -t K continuous, /(xq) = Fq, /(xat) = Fn} 

0 = {g\ 9 : [yo,yM] -t K continuous, g{yo) = Go,g{yM) = Gm} 
be metric spaces with the metrics 

dT{fi,f 2 ) = max{|/i(x) - / 2 (x)| \ x € [xo,XAr]} 

and 

dg{gi,92) = max{ \gi{x) - g 2 {x)\ \ x € [?/o,2/m] } 

respectively. Note that {F^t) and (G,dg) are complete metric spaces, hence F x Q is 
also complete with the metric d where 


d{{fi,gi), if 2 , 92 )) = m.ax{djr{fij2),dg{gi,g2)}. 


We will define a contraction mapping on F” x ^ using the coefficients , fA 

and the vertical scaling factors dA < 1 of the GIFS for a,/3 € {1,2} and 

z = 1, 2,... , iF"^. Employing the following affine invertible functions 


li : [xo,XAr] [xj_i,Xi], Ii{x) = af-x + eY 
di : [yo,yM] -t [xi-i,Xi], Ii{x) = oYkh ^ + 4-^11 
Ji : [xo,XAr] ->• [yi_i,yi], Ji(x) = af^x -h ef 
di : [yo,yM] -t [yi-i, 2 /i], Ji{x) = x 


for z = 1,..., 

for i = + 

for z = 1,..., 

for i = + 


we define the mapping 


T : Fxg F xQ 
T{f,g){x,y) = {f{x),g{y)) 
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where 


fix) 


/ 

+ dj^f{i-^{x)) + fi^ 


; X € [xi-i,Xi] for i = 1,..., 

; X € [xi-i,Xi] for i = + 1, ...,iV, 


9{y) = ' 


cfjT\y) + dff{J-\y)) + ff 


; y e [yj-i,yj] for j = 


(y) + df_i^2ig{Jj [y)) + ; y g [yj-i^yj] for j = k + 1,...,m. 

It can be easily seen that /(xq) = Fq, /(xn) = Fn, 5 ( 2 / 0 ) = Go and givn) = Gm 
using the equations ©-([HI), thus to see that T is well defined it is enough to show that / 
and g are continuous. 

/ is continuous on the intervals [xi_i,Xi] since is continuous for f = 1,..., A^. So 
we must show that / is also continuous at the points {xi, X 2 , • • •, xat-i} on which the 
function / is defined by two different ways. 

Note that for f = 1,2,..., — 1 


f{xi) = c\^I^ ^(xi) + d- ^{xi)) + 
fixi) = c}l^I~+iixi) + d}l^f{I-^}^{xi)) + flli 

are both equal Fi since I~^{xi) = xn and (Xi) = Xo by using ©. 
For i = + 1,... , N — 1 


fixi) = cf_K^iIi ^(xi) + d\l^iig{Ii ^(xi)) + 
f{xi) = c^_,_]^_^ii/j^;^(xj) + 

are both equal F) since I~^{xi) = hm and I~^i{xi) = 50 by using (l6l) . 
Finally for i = 


hxi) = cj^Ii ^(xi) + dj^fili ^(xi)) + fl^ 

Rxi) = clR_KiiI~Rixi) + dlR_j^,rg{lR^ixi)) + 

are also both equal to Fi since I~^{xi) = xn and (xi) = 2/0 by using © and © which 
shows the continuity of /. 

Similarly, one can show that g is continuous. 

To see that T is a contraction, we need to show that 


d(T(/i,5i),r(/2,52)) < r max{djr(/i,/2),dg(5i,52)} 
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for some 0 < r < 1. Let T{fi,gi) = {fi,gi) and T{f 2 ,g 2 ) = (/ 2 ,ff 2 )- Since 


max {|/i(x) -/ 2 (x)|} 

xelxo,Xj^iil 


max {|/i(x) -/ 2 (x)|} 


= max {|4^l l/i ih ^(^)) - /2 {li ^(x)) I I a; G [xi_i,Xi]} 

< M^Vj-(/i>/2 ) 

= max {|d-!^ii| \gi {l~\x)) -52 {ir\x)) | | x G [xi_i,Xi]} 

< \d^‘^\dg{gi,g2) 


where d^^ = max{|(ip| | i = 1 ,... and = max{|dP| | i = 1 ,... we get 

dT{fi,f2) < max{d^\d^^} max{djr(/i,/2), 52)}- 


Similarly, we have 

dg{ 9 i, 92 ) < max{d^\d^^} max{dj-(/i,/2), ^5(5(1,52)}, 

where d?^ = max{|d?^| | i = 1,... and d?"^ = max{|dp| | i = 1,... Using 

these results, we write 


diT{fi,gi),T{f 2 ,g 2 )) = max{djr{fi, f 2 ),dg{gi,g 2 )} < r max{djr{fi, f 2 ),dg{ 9 i, 92 )}, 


where r = maxjd^^, which is less than 1. 

From the Banach fixed point theorem, T has a unique fixed point, say {fo, 9 o), that is 
T{foj 9 o) = {fo, 9 o)- Let F and G be the graphs of /o and go respectively. Notice that for 
X G [xo,XAr] and y G [yo,yM] 


and 


fo{afx + el^) = cj^x + dj^ fo{x) + 

foia^y + eP) = c^y + dj'^goiy) + /P 


go{afx + eP) = cfx + dP/o(x) + /P 
goiafy + ef) = cfy + dfgoiy) + /f 


for 

i = 1, 


for 

i = 1,. 

..,^12 


for 


..,^21 

for 

i = 1, 

..,^22 


which imply that 


F 

G 


^11 ^12 

[J a;P(F) U [J a;P(G) 

i=l i=l 

K^l k22 


\Juf{F)\j\Juf{G). 
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Since this equations system has a solution it must he F = and G = from 

the uniqueness of the solution. Thus the graphs of the fractal interpolation functions of 
/o and go are the attractors A^ and A^ respectively. 
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